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Abstract 

Tutte conjectured in 1972 that every 4-edge connected graph has a nowhere-zero 3-flow. This 
has long been known to be equivalent to the conjecture that every 5-regular 4-edge-connected 
graph has an edge orientation in which every out-degree is either 1 or 4. We show that the 
assertion of the conjecture holds asymptotically almost surely for random 5-regular graphs. It 
follows that the conjecture holds for almost all 4-edge connected 5-regular graphs. 


1 Introduction 

A nowhere-zero 3-flow (sometimes simply called a 3-flow) in an undirected graph G = (V, E ) is an 
orientation of its edge set E together with a function / assigning a number /(e) € {1,2} to every 
e € E such that the following is satisfied. For every vertex v € V, 

X Ae) - X A e ) = °> 

eE-D+(v) eE-D - (v) 

where D + {y) is the set of all edges oriented away from v, and D~{y) is the set of all edges oriented 
towards v. 

A well known conjecture of Tutte from 1972 (see e.g. Bondy and Murty (4] (Open Problem 48) 
and Jensen and Toft [(J Section 13.3]) asserts that every 4-edge-connected graph admits a nowhere- 
zero 3-flow. This conjecture is still open. For long, it was not even known whether or not there is 
a fixed k such that every £;-edge connected graph has a nowhere-zero 3-flow (known as the weak 
3-flow conjecture of Jaeger). Weaker versions, for k > clog 2 n for n- vertex graphs, proved by Alon, 
Linial and Meshulam [T] and Lai and Zhang [7j. Recently, the weak 3-flow conjecture was settled 
by Thomassen m, who proved that every 8 -edge-connected graph admits a nowhere-zero 3-flow. 
This was subsequently improved to k = 6 by Lovasz, Thomassen, Wu, and Zhang [ 8 ]. 

It is known (see, e.g., Seymour pOj) that a graph admits a nowhere-zero 3-flow if and only 
if it has a nowhere-zero flow over Z 3 , or equivalently, an edge orientation in which the difference 
between the out-degree and the in-degree of every vertex is divisible by 3. It has also long been 
known (see [3] and [ 6 ]) that it is enough to prove the conjecture for 5-regular graphs. Thus, Tutte’s 
conjecture has the following equivalent form. 

Conjecture 1.1 (Tutte) Every A-edge-connected 5-regular graph has an edge orientation in which 
every out-degree is either 1 or 4. 

‘Supported by NSERC. 

t Supported by an ARC Australian Laureate Fellowship. 


1 



In this paper, we show that Tutte’s conjecture holds for almost all 5-regular graphs. To state this 
precisely, we say that a property of a probability space indexed by n holds a.a.s. (asymptotically 
almost surely) if the probability that it holds tends to 1 as n tends to oo (with n restricted to being 
even for odd-degree regular graphs). Using the small subgraph conditioning method of Robinson 
and Wormald [9] (see m) we show the following. 

Theorem 1.2 A random 5-regular graph G n on n vertices a.a.s. admits a nowhere-zero 3-flow, 
that is, an edge orientation in which every out-degree is either 1 or 4. 

Since it is well known that almost all 5-regular graphs are 5-edge-connected (see e.g. [T 2 ]). it 
follows that almost all 4-edge-connected 5-regular graphs have a nowhere-zero 3-flow. 

Jaeger [5j generalised Conjecture 11.11 bv conjecturing that for any integer p > 1, the edges of 
every 4 p edge-connected graph can be oriented so that the difference between the out-degree and 
the in-degree of every vertex is divisible by 2 p + 1. Similar to Tutte’s conjecture, it is known that 
the general case can be reduced to the (4 p + l)-regular case. Alon and Pralat [2] showed that the 
assertion of Jaeger’s conjecture holds for almost all (4 p+ l)-regular graphs, provided that p is large 
enough. (The lower bound for p was not optimized, but it could not be reduced to p = 1.) The 
proof used methods quite different from the present paper, involving an application of the Expander 
Mixing Lemma to an equivalent version of the conjecture. 


2 Proof of Theorem 11.2 


The pairing model for investigating properties of random regular graphs was instigated by Bol- 
lobas [3]. This consists of dn points that are arranged in n groups (called vertices) of d each, 
arranged in pairs uniformly at random. The pairs induce a multigraph in the obvious way, and we 
refer to pairs as edges. This pairing model, called V n ,d, is useful because simple graphs occur with 
equal probabilities, and the probability that it is simple for fixed d is bounded away from 0. Hence, 
to show that the random regular graph has a property a.a.s., it is enough to show that the random 
member of the multigraph corresponding to V ny d a.a.s. has the same property or is non-simple. 
(See p2J for more information on this and other claims we make about V n d-) We will work with 
orientations of (the pairs of) a pairing in V n $ in which each vertex has in-degree 1 or 4. We call such 
orientations valid. Given an orientation, vertices of in-degree 1 will be called in-vertices, and those 
of out-degree 1 out-vertices, and each point contained in an edge oriented towards an in-vertex, or 
away from an out-vertex, is called special. Moreover, a point is an in-point if the edge containing it 
is oriented towards it, and an out-point otherwise. 

Let Y = Y(n) be the number of valid orientations of a random element of V n $. It is easy to see 


that 


EU 


y 2 )5”(5n/2)! 

M(5n) 


where 

s! 

(s/ 2 )! 2 s / 2 

is the number of perfect matchings of s points. Indeed, there are ( n “ 2 ) wa Y s to select in-vertices 
(since exactly half of the vertices must be such), 5 n ways to select one special point in each vertex, 
which determines each point to be either in or out, (5n/2)! ways to pair up the points so that each 
“in” is paired with an “out”, and M(5n ) pairings in total. Using Stirling’s formula s! ~ y/2it s(s/e) s , 
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we get 


ey 


n!5 n (5n/2)! 2 2 5n /2 /25\ n/2 r 

(n/2)! 2 (5n)! V~8~J V5 ' 


(1) 


This tells us that there are plenty of valid orientations per pairing, on average. To show that 
pairings a.a.s. have at least one valid orientation, i.e. that P(Y > 0) ~ 1, a common method 
would be to estimate EY(Y — 1), show that it is asymptotic to (EY) 2 , and then apply Chebyshev’s 
inequality. As we shall see, this fails in the present case, but only just, as there is a constant 
factor discrepancy in the asymptotics. Under such circumstances, we can hope to apply the small 
subgraph conditioning method [12]. Here the first step is again to compute EY(Y — 1), and then 
also some joint moments of Y with short cycle counts (and then hope for the best). 

To estimate EY(Y — 1), consider any two orientations of the same 5-regular graph. Suppose that 
precisely k vertices are in-vertices in both orientations, and that, of these, precisely k\\ have the 
same special point in both orientations. Since the first orientation induces n/2 in-vertices, exactly 
n/2 — k vertices are in-vertices in the first orientation but out-vertices in the second one. Of these, 
suppose that k\$ (k\$ < n/2 — k) have the two special points coinciding. Similarly, there are k 
vertices that are out-vertices in both orientations; suppose that &oo of them (fcoo < k) have the two 
special points coinciding. Finally, there are n/2 — k vertices that are out- in the first and in- in the 
second orientation; suppose that &01 of them (Aoi < n/2 — k) have coinciding special points. 

It turns out that there are no additional restrictions on these parameters, other than integrality 
and non-negativity. We define 


( m 7i 'j 

I = I(n) = |(A;, k 0 o,koi, k w ,ku ) maxjkoo, k u } < k, max{k 01 ,k 10 } <--kj 


where No = N U {0}. Fix k = ( k , Aoo> &oi> k\o, ku) € I. We next calculate the number of configura¬ 
tions, i.e. pairings with two given orientations, corresponding to this vector. There are 

n! 

Aoo^oi^io^nK^ - k 0 o)'-{k ~ k n )\(n/2 -k - k m )\{n/2 -k - %))! 


ways to partition the vertices into the eight groups. There are then 


gfeoo+fcoi+feio+fell koo)+(k— fcn)+(n/2— k—koi)+(n/2—k— fcio) _ gn _ ^n— fcoo —fcoi — fcio — 


ways to assign special points in the two orientations. Next, we need to pair (in,in)-points with 
(out,out)-points (where the first “in” refers to the first orientation, and so on), and (in,out)- with 
(out,in)-points. The number of (in,in)-points is equal to 


ku + 4Aoo + 3 (k - /eqo) + ( n/2 - k - k 0 1 ) + (n/2 - k — k w ) = n + k + k 00 + k u ~ k m - k w , 


and the same applies for (out,out). These two sets must be paired with each other. Half of the 
remaining points, or 

bn/2 - (n + k + k 00 + ku ~ km - &i 0 ), 
will be in-out, and an equal number will be out-in. Hence, there are 


(n + k + k 00 + k u ~ k 0 1 - Ai 0 )!(3n/2 - {k + fc 0 o + - ^oi - &io))! 

ways to legally pair the points. The number of configurations is the product of the above factorials. 
To obtain the expected number of pairs of orientations, we must divide by the number M(5n) of 
pairings. Putting z = z(k) = k/n and applying Stirling’s formula again, we can write 

EY(Y - 1 ) = Y1 r ( z )5 , ( z ) ex P («/(z)), (2) 

k el 
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where the various factors are defined as follows. The function r, which is the error factor in the 
applications of Stirling’s formula, has the property that r = 0 ( 1 ) for all z, and r ~ 1 if z is bounded 
away from the boundary of 

J : = j(2,z 00 ,201,210,211) € Mq : 2 < max-j>oo, 2 n} < 2, max{z 0 i,2i 0 } < ^ - z\ 

where Mo is the set of non-negative reals. With 6 = 2 + 1 + 200 — 201 — zio + 211 and h(x) = x log x, 

1 / 6(5-26) \ 1/2 

g v / 32 ( 7 T 77,) 5 / 2 v200201210211(2; - z 0 o)(2 - 2n)(l - 22 - 2 zi 0 )(l - 2 z - 2 z 0 l)J 

from the polynomial factors in Stirling’s formula, and 

/ = ( 9/4 - Zoo - 201 - 210 - 211 ) log4 + log5 - 6(5) + 6(5/2) + h(b) + 6(5/2 - 6) - h(z 0 0 ) 

-h(zoi) - 6,(210) - 6(211) - 6(2 - 2 00 ) - 6(2 - 2n) - 6(1/2 - 2 - zqi) - 6(1/2 - 2 - 210) 


from the rest. 

Note that we can extend the definition of / continuously to the boundary of J by defining 
xlogx = 0 at x = 0. Then / achieves its maximum on J. Our next goal is to show that z = 
(1/4,1/20,1/20,1/20,1/20) is the unique global maximum point of / on J, since we can easily 
argue then that points far away from z give negligible contribution to (|2|). 

We first investigate stationary points in the interior of J. An algebraic manipulation package, 
such as the Maple we used, makes this easy. We find 

df = 1q (2 - 2 00 )(2 + 1 + 2 00 - 2qi - 2iq + 2n) 

dzoo ^ 22oo(3 — 2 2 — 2zoo + 2201 + 2z\o — 2zu) 

Setting this equal to 0 gives Poo = 0 where 

Poo = (z - 2 00 )(2 + 1 + Zoo - Z 01 - 210 + 211 ) - 2200(3 - 22 - 2200 + 2201 + 2210 - 2211 ). 

Defining P 01 etc. similarly, and P from df/dz, we obtain five polynomials such that any local 
maximum must be a common zero of all five polynomials. Write Y, x) for the resultant of two 

polynomials X and Y with respect to x. When X = Y = 0, it is necessary that &(X,Y,x)= 0. We 
find that 


P 6 := ^(P 00 , P01,210) = 52 - 52 00 - IO2 2 + IO2200 - IO2012 - 1502 0 i 2 00 , 

P 6 - ^(Poo,Plo,2oi) = 10(-2io + 201))(2 + 15 z 00 ). 

On the interior of J, we have 2 > 0, so we may assume that 201 = 210 • Also, 

^(P 6 ,^(Poi,P n ,2io),2oi) = -8002(200 - 2n)(-l + 2 2). 

On the interior 2 < 1/2, so we may assume that 200 = Z\\. 

For any polynomial X , let X* be the result of setting 2 n = 200 and 210 = 201 . We find 

P 7 := ,^(P 0 * 0 , P 0 *,, 2 01 ) = -6O2 3 - 480 z 2 2 00 - I2O2200 - 150022 ^ + 2040 2^ - 18002 ^ 0 , 

^ , (P 7 ,^ , (P*,Poo, 2 oi), 2 oo) = 573308928 x 10 11 x 2 8 (42 - 1)(130682 2 - 65342 - 109)(1 - 2 zf. 
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Neither root of the quadratic factor is in [0,1/2], so at any interior stationary point we must have 
z = 1/4. Substituting z = 1/4 into P 7 = 0 gives 

(20z 0 o - 1)(96 zq 0 - 84z 00 - 1) = 0. 

Again, the roots of the quadratic factor are out of range, and hence zoo = 1/20 at any interior 
stationary point. Thus zn = 1/20 also. Substituting the known values into Pq X gives 

^(20z 0 i-l)(24zoi-23) 

and hence the unique stationary point of / in the interior of J is z. 

We turn next to the boundary of J. First consider any point on the boundary at which 0 < 
z < 1/2. Then zqi + zio < 1 — 2z and zoo + 211 < 1 < 3/2 — z. Hence, for zoo tending towards 0 
or z, df/dzoo is dominated by the terms — logzoo and +log(z — zoo), which tend to +00 and —00 
respectively in the two cases. This shows that for 0 < z < 1/2, the boundary points where zoo = 0 
or zoo = 2 cannot contain a global maximum of / on J. A similar observation applies also to show 
that zn, zio and zoi cannot be at the extreme ends of their ranges. Hence, no such boundary point 
is a maximum of /. 

We are left with considering the boundary points where z = 0 or z = 1/2. 

Case 1. z = 0. 

Membership of z in J then forces zoo = 211 = 0. Substituting z = zoo = zn = 0 into / gives a 
function /(zoi,zio) with domain [0,1/2] 2 . Setting its partial derivatives to 0 give two polynomials 
equal to 0. Their difference is 

-2(6z 0 i - 11 + 6zi 0 )(-zio + zoi). 

Since zoi + 210 < 1, this is only 0 at zoi = zio- We find df(z\o,z\o)/dz\o = (3 + 4 zoi) 2 / 256 zq 1; 
which is never 0 . So no interior point can be maximum for /. This leaves the boundary of its 
domain, where zoi or zio is 0 or 1/2. If zoi = 0 , then dfdzoi is large and positive, so there is 
no maximum there. By symmetry, it is a similar conclusion if zio = 0 . On the other hand, along 
the boundary zoi = 1/2, d/(l/2, zio)/dzio = log((2 + zio)/ 4 zio), which is always positive, and 
similarly for zio = 1/2. So the only possible local maximum on the boundary in this case is at 
( 0 , 0 , 1/2,1/2, 0 ). Here / = log(5/8), whereas /(z) = log(25/8). 

Case 2. z = 1/2. 

Membership of z in J then forces zoi = Zio = 0. It turns out that substituting z = 1/2 and 
zoi = z 10 = 0 into / produces exactly the function /(zoo, 2 n), with the same domain. Hence Case 1 
shows that the only local maximum of / here is log(5/8). 

In conclusion, / has (at most) three local maxima, at z, (0,0,1/2,1/2, 0) and (0,1/2,0,0,1/2), 
and the first of these is the global maximum on J. 

The next part of the proof consists of a routine computation and argument. Set z = 1/4 + y and 
Zij = 1/20 + yij for each of the other four variables. Then expand /(z) about the global maximum 
point z, to obtain 

/ = log(25/8) + y T By + 0(x 3 ) (4) 

where y = (y, y 00 , yoi, Vw, Un), y T denotes the transpose of y, 


-92 

33 

-33 

-33 

33 

33 

-117 

-8 

-8 

8 

-33 

-8 

-117 

8 

-8 

-33 

-8 

8 
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-8 

33 

8 

-8 

-8 
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and x = x(z) = ||y||, with || • || denoting the L 2 norm (say). The error term in this expansion is 
valid by Taylor’s theorem provided that x = o(l). 

A standard argument (with details given below) now allows us to estimate the summation in ([2]) . 
and we obtain 


E Y(Y - 1) ~ 


25^ n 5'(z)(7m) 5//2 
8 ) y/\ det B\ 


(5) 


We find that g{ z) = (5 5 /2)(7rn) 5//2 and det B = —3• 5 6 - 7/4, so P) gives ET(T —1) ~ -^=(25/8) n . 
Combining this with (JT|), we have 

EY 2 5 

Wy ~ 7Ti' (6) 

We now justify (f5|) . Let Jo := {z : x = o(n -2 / 5 )}. For z G Jo, we have r(z)g(z) ~ g( z) and 
x 3 = o(n _6//5 ). Thus 


Y r(z) 5 (z)exp (n/(z)) ~ fif(z) Y e nyTBy . 

k:k/n£Jo k:k/raGJo 


The eigenvalues of B are (—37 ± \/697)/4, each with multiplicity 1, and —25/2 with multiplicity 3. 
These are all less than —2.6. Hence B is negative definite, and we have 

Y e nyT ' By ~ n 5 f e nyT By d n y. 
k:k/nSJ 0 ^ J ° 


Here the factor n 5 arises from the change of variable from k to z = k/n, and the replacement of the 
summation by the integral can be justified by various elementary means. (One way is to rescale to 
w = v/ny and observe that the summations for various n amount to different Riemann sums for 
the integral of the fixed function e w . The regions of integration are different for different n but 
the tails of the summation and the integral outside the regions are easily seen to be negligible.) A 
simple argument, for example bounding the value of the integrand in terms of ||y||, shows that the 
latter integral can be extended to all of M 5 with no significant change, that is 



e ny T By d n y 





e nyT By (jriy = 


7T 5 / 2 

n 5 / 2 -^/| det B\ 


Combining the last few conclusions gives 


Y r ’( z )^( z ) ex p ( n /( z )) ~ 

k:k/n£ Jo 


25^” (/(z^tto) 5 / 2 

8 J yj\ det B\ 


(7) 


Recalling that B is negative definite, and we conclude from (J4j) that on the boundary of Jo, 
the value of / is /(z) — H(n” 4//5 ). Since / is a fixed function independent of n, and z is its global 
maximum, 

max f(z) = /(z) — n(n _4/5 ) = log(25/8) — H(n -4 ^ 5 ). 
zGJ\Jo 

Considering that r and g are polynomially bounded, it follows that each term in the summation 
in (|2]) (where / is multiplied by n) for which k/n € J\ Jo is (25/8) n e _J2( ' nl/,J ' ) . There are a polynomial 
number of such terms, so their sum is likewise (2h/&) n e~^ nl/ Thus the terms in (0) dominate 
the summation in Q, and we have (J5]) as claimed. 
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The second main computation in the small subgraph conditioning method involves the random 
variable X k (k > 1) which is the number of cycles of length k in V n $. It is known that for each 
k > 1, X \, X 2 ,... ,X k are asymptotically independent Poisson random variables with 


^X k ~^ — 


4 fc 

2 k' 


We are required to show, for each k > 1, that there is a constant /i^ such that 

E (YX k ) 

E y -> Vk 

and, more generally, such that the joint factorial moments satisfy 

‘' [Xk]j k ) 

EY 11 fi 


( 8 ) 

(9) 


( 10 ) 


for any fixed j\ ,..., j k - (Here [ x\k is the falling factorial.) We will derive a value of /i^ satisfying Q 
and will observe that essentially the same argument generalises easily to give m- 

To evaluate E (YXk), we find the number of triples ( P,C,0 ) where P is a pairing, C a fc-cycle 
of P and O an orientation of P, and then divide by \P n $\ = M(5n). 

The number of ways to choose the pairs of (i.e. inducing the edges of) the cycle is 20 k [n] k /2k. 
The easiest way to see this is to choose the vertices of the cycle in one direction, starting at a 
canonical one ([ra]fc), and then an ordered pair of points in each vertex (5 fc 4 fc ), and divide by the 
multiplicity of counting due to the canonical start and direction (2k). 

We will count the triples ( P , C, O ) which have i vertices on C with in-degree 2 in C. These must 
be out-vertices in the orientation O. The number of vertices with out-degree 2 in C must also equal 
i, and these must be in-vertices. The remaining k — 2i cycle vertices have in-degree and out-degree 
in C both equal to 1. These can be either in- or out-vertices. We now select the rest of the oriented 
pairing. There are (^-J wa Y s 1° select the remaining in- and out-vertices, and 3 2 * ways to choose 
the special points of the vertices of C. (Note that the latter only needs to be done for vertices of 
in-degree 0 or 2 in C\ vertices of in-degree 1 in the cycle have their special point already determined 
by the orientation of the edges of the cycle.) Finally, there are 5 n ~ k ways to choose the special 
points of vertices outside C, and then (5n/2 — k)\ ways to pair up the points of appropriate types. 
Let a,i denote the number of orientations of C with i vertices of in-degree 2. Dividing the number 
of triples by M(5n)EY, and summing over all i, we obtain 


E (YX k ) 
E Y 




£ 

0<i<k/2 


20 k (n} k {^ 2 2 ! l )3^- k (5n/2 - k)\ 
2fcy 2 )5-(5n/2)! 



This gives Q with 



To find a*, one can select the 2 i vertices of C that are to have out-degree 0 or 2 in C, and after 
this there are exactly two ways to orient C. Hence dj = 2(,^ : ), and this is the coefficient of x 2l in 
q(x) := 2 (1 + x) k . It follows that 


2 ai 

0<i<k/2 


2 i 


+ <?(- 
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and thus 


^ = 2fc( 4fc + ( " 4/5)fc )- 

As mentioned above, the derivation of (j9|) is a straightforward generalisation of this. One starts 
with a set of cycles instead of one cycle, and following the same argument, the effects of those cycles 
lead to multiplicative factors in the counting that turn out to be asymptotically independent. We 
omit the details. 

The final step of small subgraph conditioning requires us to compute (see ([ 8 ])) 




Rk_ 

A k 



k 


and then, using — log(l — x) = J2k>i x k /k, 

ex p ( g Wl) = exp (i g i (|) ) = exp (-1 log (‘ - 1 ) ) = -|=. 

The fact that this is coincides with the right hand side of ([ 6 ]) implies, by m Theorem 4.1], that 
P(Y > 0) ~ 1. That is a random multigraph in V n .d a.a.s. has a nowhere-zero flow over Z 3 . 
Theorem Ol now follows, in view of the comments at the start of this section. 


References 

[1] N. Alon, N. Linial, and R. Meshulam, Additive bases of vector spaces over prime fields, J. Com¬ 
binatorial Theory, Ser. A 57 (1991), 203-210. 

[2] N. Alon and P. Pralat, Modular orientations of random and quasi-random regular graphs, 
Combinatorics, Probability and Computing 20 (2011), 321-329. 

[3] B. Bollobas, A probabilistic proof of an asymptotic formula for the number of labelled regular 
graphs, European Journal of Combinatorics 1 (1980), 311-316. 

[4] J.A. Bondy and U.S.R. Murty, Graph Theory with Applications. Macmillan, London, 1976. 

[5] F. Jaeger, Nowhere-zero flow problems, in: L. Beineke, et al. (Eds.), Selected Topics in Graph 
Theory, vol. 3, Academic Press, London, New York, 1988, pp. 91-95. 

[ 6 ] T.R. Jensen and B. Toft, Graph Coloring Problems. Wiley-Intersci. Ser. Discrete Math. Optirn., 
1995. 

[7] H-J. Lai and C.Q. Zhang, Nowhere-zero 3-flows of highly connected graphs, Discrete Math. 
110 (1992) 179-183. 

[ 8 ] L.M. Lovasz, C. Thomassen, Y. Wu, and C.Q. Zhang, Nowhere-zero 3-flows and modulo k- 
orientations, Journal of Combinatorial Theory, Series B 103 (2013) 587-598. 

[9] R.W. Robinson and N.C. Wormald, Almost all cubic graphs are hamiltonian, Random Struc¬ 
tures Algorithms 3 (1992), 117-125. 

[10] P. D. Seymour, Nowhere-zero flows, in “Handbook of Combinatorics,” 299, North-Holland, 
Amsterdam, 1995. 

8 



[11] C. Thomassen, The weak 3-flow conjecture and the weak circular flow conjecture, Journal of 
Combinatorial Theory, Series B 102 (2012) 521-529. 

[12] N.C. Wormald, Models of random regular graphs, Surveys in Combinatorics, 1999, J.D. Lamb 
and D.A. Preece, eds. London Mathematical Society Lecture Note Series, vol 276, pp. 239-298, 
Cambridge University Press, Cambridge, 1999. 


9 


